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TECHNICAL NOTE 3888 

SIMPLIFIED METHOD FOR ESTIMATING COMPRESSIBLE LAMINAR HEAT 
TRANSFER WITH PRESSURE GRADIENT 
By Eli Reshotko 


SUMMARY 

In the present report an approximation is made in the method of Tech- 
nical Note 3326, -which simplifies the calculation of heat transfer. Good 
agreement with the method of Technical Note 3326 is expected for isother- 
mal surfaces with adverse and small favorable pressure gradients regard- 
less of surface temperature and also for flow with large favorable pres- 
sure gradient over highly cooled surfaces . 


INTRODUCTION 

A major factor influencing the design of vehicles for supersonic and 
hypersonic flight is aerodynamic heating. While for simple shapes such 
as the flat plate, the cone, and stagnation points the amount of laminar 
heat transfer may be estimated using simple formulas, such is not general- 
ly the case for a calculation with arbitrary pressure gradient. 

An integral method for calculating compressible laminar boundary 
layers with heat transfer and pressure gradient is given in reference 1. 
This method is an adaptation of the correlation technique of Thwaites 
(ref. 2) to the compressible boundary layer using the exact solutions of 
reference 3 . The present report points out an approximation which sim- 
plifies the calculation of laminar heat transfer. Examples are presented 
to illustrate the approximate technique, and numerical comparison is made 
with the complete method of reference 1. Some approximate techniques for 
slender bodies introduced by Jack and Diaconis (ref. 4) are described in 
conjunction with the examples . 


SYMBOLS 

A coefficient in eqs. (2) and (3) 

B exponent of Mach number in eqs. (2) and (3) (fig. l) 
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C f Re w 

Nu 


D 

h 

K 

k 

L 


M 

Nu 

n 


P 

R 

Re 


local skin-friction coefficient, 


yi 


pw^i 


Reynolds analogy parameter (ref. l) 

P - Poo 

pressure coefficient, ^ 

STPcoMl 

specific beat at constant pressure 

throat diameter of rocket nozzle (example iv) 
heat-transfer coefficient 

a f jrr-rr < K ' 4 ft>r r - 1 -*) 

thermal conductivity 
arbitrary length 


dimensionless shear parameter (ref. l) 


9 t v ( du\ 

' u e t e W/v 


Mach number 


lix 

Nusselt number, r— 

% 


correlation number 




dimensionless pressure gradient 


P C T 


i ^/ rM * 

' P e ax / rM e 


Pr Prandtl number, -rr— (evaluated at wall temperature) 


pressure 

radius of axially symmetric body 
Pw u e x 


w 


Reynolds number, 




V 


t 


» 
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r / S v 

t 

U 1 


nflt w [ay(t/t e )] w 



temperature 
external velocity 


x 

y 

a 

r 

e 

p 

V 

p 


coordinate along surface 
coordinate normal to surface 

exponent of Prandtl number for heat -transfer modification 

ratio of specific heats 

houndary-layer momentum thickness 

absolute viscosity 

kinematic viscosity 

density 



Subscripts: 
cyl cylinder 

e local flow outside boundary layer (external) 

w wall or surface value 

“ free-stream static value 

0 free-stream stagnation value 


CALCULATION PROCEDURE 

The calculation of laminar heat transfer over isothermal surfaces by 
the method of reference 1 is first summarized and then the approximation 
introduced . 

The heat -transfer parameter Nu/^ Re v for Prandtl number 1 from 
equations (36) and (37) of reference 1 is 
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The factor 




which is directly proportional to the square of mom- 


entum thickness, is obtained in reference 1 by quadrature from the momen 
turn integral equation. For two-dimensional flow 



( 2 ) 


while for axially symmetric flow 



Reference 1 suggests that for moderate pressure gradients the coefficient 
A be set equal to 0.44 and the coefficient B be chosen from a figure 
which is presented herein as figure 1. If an expected favorable pressure 
gradient is consistently that corresponding to stagnation-point flow or 
stronger, the coefficients A and B should be evaluated from the ap- 
pendix. Such large pressure gradients might be encountered, for example, 
on the nose portions of blunt bodies or In choked nozzles . 


Thus, all quantities in equation (l) are determined directly by the 
Mach number and static-temperature distributions with the exception of 
the shear or skin-friction parameter l and the Reynolds analogy param- 
/Cf-Rejv 

eter | -= ) . These parameters are presented in reference 1 as func- 

\ Nu /Pr=*l 

tions of the correlation number n. In order to complete the calculation 
of heat transfer according to reference 1, the parameter n must be eval- 


uated from the quantity 



which requires differentiation of the 
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\ Ku / Pr=3 


pressure distribution; the parameters 2 and ^ — / are then ob- 

*Pr=l 

tained and substituted into equation (l). However, as can be seen from 
figure 2, the ratio ^ 




which is the desired combination of 


Pr=l 

parameters, is relatively insensitive to 
Hence, it is suggested that 


m 


n for many cases of interest, 
be taken as 0.11 for all values 


'Pr=l 

of n. 1 This assumption gives an error of less than 10 percent except 
for flows with large favorable pressure gradients or very near separation. 
For large favorable pressure gradients, a further suggestion is made in 

the appendix with regard to the quantity ^ 


/Wi ■ 

/Pr=l 


The expression for heat transfer upon setting 


0.11 is 


( C f* e w\ 

\ Nu /Pr*l 


equal to 


where 


Pr 0, accounts for the effect of Prandtl number on heat transfer. Values 
of a suggested by Tifford and Chu (ref. 6) are, for small pressure 
gradients, l/3; for large adverse pressure gradients, l/4; and for ex- 
treme favorable pressure gradients, l/2. 

It should be noted and emphasized here, that if skin-friction in- 
formation is desired, the parameter 2 must be evaluated, and, therefore, 
the complete method of reference 1 must be used. 


■bi 


(v^) 


0.22 Pr a 


x 

L 


1 p, $ 


( 4 ) 


is obtained from either equation (2) or (3) . The quantity 


^This same approximation can be obtained from the almost linear vari- 
ation of a heat-transfer parameter r/s^ with correlation number n 

shown in fig. 6 of Crocco and Cohen (ref. 5). The present approximation 
corresponds to r/n S^. = -0.22. 
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EXAMPLES 

Three examples are given. The first two are a cone cylinder and a 
parabolic -nose cylinder for which Jack and Diaconis (ref. 4) recently 
presented experimental information. The numerical calculations for these 
two examples are for a free-stream Mach number of 3.12, ambient static 
temperature of 180° R, and for the surface cooled to ambient static tem- 
perature. The third example is that of the rocket nozzle calculated in 
reference 1. 


I. Cone Cylinder 


Because the cone cylinder is composed of two simple shapes, it lends 
itself to a somewhat generalized treatment. The length of the cone sur- 
face will be designated L as shown in the following sketch: 



Over the conical portion of the body the surface Mach number is constant 
so that from equation (3) 




and from equations (4) and (5) 


( \ = 0.575 Fr 1 / 3 

\ / V Re w/cone 


(5) 


( 6 ) 


in which the coefficient 0.575 is identically that indicated by the 
exact boundary-layer solution for a Prandtl number of 1. 
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On the cylindrical portion of the "body from equation (3), 



and the heat-transfer parameter Fu//^Re w is obtained by substituting 
equation (7) into equation (4). 

The heat-transfer rate on the cylinder side of the cone-cylinder 
junction can be readily evaluated. From equations (4) anri ( 7 ), 



Equations (5) to (8) are used to calculate the heat transfer to the 

1 ° 

-included-angle cone cylinder reported in reference 4. The pressure 

distribution for this cone cylinder, shown in figure 3, was obtained using 
the theo ry of Van Dyke (ref. 7). An examination of the calculated values 
of Nu/ 'V^ e - W - shows that the presented approximate technique gives results 

within 2 percent of those using the complete method of reference 1. The 
abrupt change in the heat-transfer parameter at the cone-cylinder j unc tion 
is due to the abrupt break in the pressure and Mach number distribution at 
this station. 


Jack and Diaconis (ref. 4) have proposed that the heat tr ans fer on 
the cylindrical portion of the body be estimated by a formula which 
matches the cylinder value at the cone -cylinder junction and the flat- 
plate value at infinity. This relation takes the form (for x > L) 
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where 

^ _ \/>/Re w /cyi at junction 

(JO?—) 

' cone 

as evaluated from equation (8) . The use of equation (9) avoids the neces- 
sity of performing the integration of equation (3) over the cylindrical 
portion of the "body., and it is seen from figure 3 that equation (9) agrees 
closely with the othef estimations of heat transfer on the cylindrical 
portion of the "body. All three curves in figure 3 represent the data of 
reference 4 to within 15 percent. 


II. Parabolic Nose Cylinder 

A sketch of the parabolic nose cylinder together with its pressure 
distribution as calculated using the theory of reference 7 i s sho wn in 
figure 4. The distribution of heat-transfer parameter Nu//^Re w as cal- 
culated by the present approximate technique using equations (3) and (4) 
is seen to be within 3 percent of that obtained using the complete method 
of reference 1. 

A further simplification for- this slender parabolic-nosed body is 
proposed in reference 4. If it is assumed that the Mach number over the 
body is fairly constant and that the heat -transfer variation is due prim- 
arily to the geometric variation of the body from conical to cylindrical 
shape, then from equations (l) and (3) 



The curve of equation (10) is seen from figure 4 to be in good agree- 
ment with those of the complete (ref. l) and simplified techniques. All 
three curves represent the experimental data of reference 4 within 12 
percent . 


III. Rocket Nozzle 

The third example is that of the rocket nozzle (fig. 5) also used as v 

an example in reference 1. It is an example involving large favorable 
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pressure gradients with heat transfer. The combustion-chamber stagnation 
pressure is assumed to "be 500 pounds per square inch absolute, the stagna- 
tion temperature is taken as 4000° E, the Prandtl number is assumed to be 
0.78, and the ratio of specific heats is 1.3. The nozzle wall is assumed 
cooled to a uniform temperature of 800° E, which corresponds to 

t w/ t 0 = °' 2 * 

The solid line in figure 5 is the variation of the heat-transfer 
parameter Nu/^/Ee^ from the complete method of reference 1. The sim- 
plified technique was applied in two ways: first, as proposed in the 

text of this report and second, as proposed for large favorable pressure 
gradients in the appendix. The results are shown in figure 5. While 
both applications of the simplified technique result in an overestimation 
of the heating lead, the technique for large pressure gradient matches the 
result of reference 1 at the stagnation point and overestimates the heat- 
transfer rate at the throat only by about 13 percent. 


CONCLUDING EEMAEKS 

An approximation is made in the method of Technical Note 3326 which 
simplifies the calculation of heat transfer. 

The simplified method is expected to be adequate for isothermal sur- 
faces with adverse and small favorable pressure gradients regardless of 
surface temperature level and also for flow with large favorable pressure 
gradient over cooled surfaces . 

Good agreement is obtained in these ranges between the simplified 
technique and the complete method of Technical Note 3326 for the illus- 
trative examples. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, September 13, 1956 
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APPENDIX - CALCULATION PROCEDURE FOR LARGE FAVORABLE PRESSURE GRADIENTS 


The simplified, technique as proposed in the text of this report is 
set up in such a way as to match the exact "boundary-layer solution for 
zero pressure gradient. For large favorable pressure gradients (corres- 
ponding to stagnation-point flow or stronger), the use of the values of 

A, B, and -jtts \ proposed in the text may greatly overestimate the 

/ L f Re v \ 

\ Wu /Pr=*l 

amount of aerodynamic heating. 


The values of A and B in the following table correspond to lines 
tangent to the curves of figure 3 of reference 1 at the plane stagnation 

l 


point. Similarly, the values of 


nation point flow. 


C 


Cf Re w\ 


and a are for plane stag- 


Nu 


'Pr=>l 


t w 

A 

B 

l 

a 

t 0 

/ C f Re w\ 

\ Nu / Pr=l 

0 

0.380 

2.35 

0.102 

0.4 

.2 

.372 

2.53 

.099 

.4 

.6 

.362 

3.20 

.093 

.4 

1.0 

.410 

4.94 

.083 

.4 


The simplified technique presented in this report is not recommended 
for use with large favorable pressure gradients over heated surfaces. 

While the numbers presented in this appendix are for matching at 
plane stagnation-point flow, a similar match could be made at any desired 
condition using the information in figure 3 of reference 1 and in figure 
2 of thiB report. - 
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